Tetthetsfunksjonalteori

Dette notatet er en kort introduksjon til tetthetsfunksjonalteori, eller "density-functional theory" (DFT) på engelsk. Tetthetsfunksjonalteori er kvantekjemiens ubestridte arbeidshest. De aller fleste beregninger i kjemi gjort i dag er med DFT. Grunnen er at det er den eneste metoden som kan brukes med dagens datamaskiner og samtidig regne ut noenlunde nøyaktige egenskaper for molekyler av interessant størrelse. DFT er mer nøyaktig enn Hartree-Fock, mens kostnadene er omtrent de samme.

NB: DFT er pensum. Dette notatet er ment som en orientering og generell beskrivelse, og det ikke meningen at du skal kunne gjenta utledningene her. De er uansett ikke helt presise, og mange tekniske steg er forbigått. Det viktigste er at DFT kan ses på som en korreksjon til Hartree--Fock, og at Kohn--Sham-likningene har samme struktur som Hartree--Fock. Det er viktig å opparbeide seg en viss forståelse for hvordan Kohn-Sham-likningene er bygget opp, og hva det teoretiske grunnlaget for DFT er.

Walther Kohn

I 1999 fikk Walther Kohn nobelprisen for oppfinnelsen av DFT. (Prisen ble delt med John Pople som fikk pris for sitt bidrag til kvantekjemi.) Han var en Østerriksk-Amerikansk jødisk fysiker, som ble reddet fra Østerrike til England av de allierte rett etter at Hitler hadde annektert landet, i 1938-1939. Han ble fraktet til Toronto i Canada. Han fikk ikke lov å studere kjemi siden han var tysk, og studerte derfor fysikk og matematikk. Likevel fikk han i 1998 nobelprisen i kjemi for sitt arbeid med DFT! Han tok PhD ved Harvard i USA, og ble professor ved University of California at Santa Barbara.

walther-kohn.png

Hohenberg-Kohn-teoremet og dets konsekvenser

Schrödingerlikninga, en påminnelse

For et gitt molekylært system har vi et ytre kjernepotensial

Her er kjerneposisjon til atom nummer , og er kjerneladningen i enheter av .

Vi søker grunntilstanden Det vil si at vi har en løsning av Schrødingerlikningen med minst energi :

Ekvivalent ønsker vi å minimere forventningsverdien til energien over alle normerte bølgefunksjoner. (Variasjonsprinsippet/metoden.)

Vi bemerker at bølgefunksjonen er en svært komplisert funksjon og at vi er motivert til å vinne en enklere beskrivelse.

Elektrontettheten

Hohenberg-Kohn-teoremet er en bemerkelsesverdig observasjon som omhandler elektrontettheten, en av de viktigste observablene. Vi definerer tettheten assosiert med en -elektron-bølgefunksjon: Tettheten er ganger sannsynligheten for å finne minst ett elektron i posisjonen . Formelen er

der vi for enkelhets skyld har sett bort fra at elektronet har spinn. Å ta med spinn gir bare en mer komplisert formel med mindre intuisjon, uten å endre resultatet. Dersom vi integrerer over hele rommet, så ser vi at vi får

som stemmer med at representerer en tetthet. Elektrontettheten er hva vi mener når vi snakker om "elektronskyen".

Første Hohenberg-Kohn-teorem

Det første HK-teoremet sier kort sagt at dersom er tettheten som kommer fra en grunntilstand for et molekyl med kjernepotensial , da er bestemt opp til en konstant , dvs., alle andre potensialer (ikke bare kjernepotensialer) som er grunntilstand for, er på formen . Altså, dersom vi flytter på kjernene, selv bare ørlittegrann, så vil også tettheten endre seg.

Vi observerer at og har den samme grunntilstandsbølgefunksjon, og at det er bare ett valg av som gir et ekte kjernepotensial, nemlig .

Vi kan nå gjøre følgende observasjon: Gitt en tetthetsfunksjon , så bestemmer denne funksjonen kjernepotensialet (opp til en konstant som vi kan sette til null). Men dersom vi kjenner kjernepotensialet, kan vi - i prinsippet - løse Schrödingerlikninga og finne grunntilstanden til dette potensialet. Men da kan vi regne ut alle forventningsverdier, og dermed all fysikk!

Med andre ord: Det er nok å vite tettheten for å bestemme løsningen av Schrödingerlikninga!

Andre Hohenberg-Kohn-teorem

Det andre HK-teoremet sier at det finnes en funksjonal , kalt den universelle funksjonalen, slik at grunntilstandsenergien til kan skrives

En funksjonal er en funksjon av en funksjon. Tenk på en funksjonal som en funksjon som tar uendelig mange varable som innputt, og produserer ett eller flere tall som utputt. Vi har understreket i notasjonen på venstre side at grunntilstandsenergien er en funksjonal av potensialet.

Likningen kalles ofte "Hohenberg-Kohn-variasjonsprinsippet." Funksjonalen er "universell" i den forstand at den er gyldig for alle molekyler. Merk at vi har fiksert antall elektroner i denne sammenhengen, men prinsippet kan utvides til alle .

Levy-Lieb constrained-search-argumentet

Eksistensen av en universell funksjonal kan virke overraskende, men vi kan bruke et enklere argument som ikke tyr til HK-teoremene for å vise at den må finnes.

Vi har dette uttrykket for forventningsverdien til Hamiltonoperatororen:

Her er tettheten assosiert med , og vi understreker her avhengigheten med subskriptet.

Nå observerer vi at vi skal minimere denne forventningsverdien over alle normerte bølgefunksjoner.

Det er mange bølgefunksjoner som git opphav til samme tetthet. Men alle disse bølgefunkjonene gir samme forventningsverdi av ! Så vi deler opp minimeringen over i (1) minimering over og (2) minimering over alle som har tetthet . Dette siste skrives . Da får vi:

Den siste likninga definerer en Levy-Lieb-versjonen av den universelle funksjonalen. Den som nå er klart er at den representerer den kinetiske energien og elektronvekselvirkningene til faktisk alle bølgefunksjoner slik at .

Sammenliknet med (fra HK2) er eksplisitt definert. De to funksjonalene er matematisk ikke helt identiske, men gir opphav til samme grunntilstandsenergi.

De matematiske detaljene angår oss ikke. Vi vil heretter snakke om "den universelle funksjonalen ", og ikke bry oss om hvilket prinsipp den kommer fra. Den eksisterer!

Den universelle funksjonalen er en ekstremt komplisert funksjon som avhenger av i alle punkter i rommet på en ekstremt innfløkt måte. Den er så innfløkt at vi ikke vet helt hvordan den ser ut. Det er rett og slett et mysterium! Matematikken gir oss i essensen stort sett bare at denne funksjonen eksisterer, og noen svært få hint om hvordan den eller ser ut. Men vi må gjøre kvalifiserte gjetninger og approksimasjoner for å komme noen vei videre.

Oppsummering av fundamentet for DFT

Den veldig kompliserte -elektron-bølgefunksjonen er byttet ut med en mye enklere tetthet, og at Hamiltonoperatorens komponenter er byttet ut med den ukjente (og mye mer kompliserte) funksjonalen . Grunntilstanden er gitt ved et variasjonsprinsipp.

Kohn-Sham-teori

Det er ikke meningen at du skal forstå denne "utledningen". Den er til orientering og målet er å gi en viss ide om hvor Kohn-Sham-likninga kommer fra, og hva de ulike symbolene betyr.

Et vanskelig problem med DFT er at er så utrolig komplisert og vanskelig å modellere. Særlig to aspekter ved er vanskelige å kontrollere: Det er vanskelig å konstruere kinetisk energi fra bare tettheten alene, og det er vanskelig å modellere effekten av Pauli-prinsippet, "exchange", samt vekselvirkningene mellom elektronene, kalt "correlation" (korrelasjon på norsk, men vi bruker engelsk term her).

I Hartree-Fock er begge aspekter lette å modellere - kinetisk energi er jo bare en operator, og pauliprinsippet gir "exchange-potensialet" . Potensialet kalles "Hartree-potensialet" og er faktisk ikke vanskelig å modellere. Det er eksplisitt uttrykt ved tettheten ! Hartree-potensialet er det elektrostatiske potensialet generert av ladningstettheten , og velkjent i elektrodynamikken:

Den totale elektrostatiske energien er

Kohn og Sham hadde den briljante ide at man kunne gjenta HK-teoremene for et ikkevekselvirkende system av elektroner, og betrakte et slikt - fiktivt - system i et annet ytre potensial . Dette potensialet tenkte de seg kunne velges slik at grunntilstanden for det fiktive systemet reproduserte den eksakte vekselvirkende tettheten. Hvorfor ikke-vekselvirkende? Jo, for da stryker vi bare i Hamiltonoperatoren, som gjør Schrödingerlikninen lett å løse ved separasjon av variable, som gir

Den eksakte ikke-vekselvirkende -elektron-tilstanden er en Slater-determinant, akkurat som i Hartree-Fock-teori.

Ved hjelp av HK-teoremene eksisterer det nå en universell funksjonal for ikkevekselvirkende systemer, kalt , slik at den ikke-vekselvirkende grunntilstandsenergien er

Men samtidig vet vi fra separasjon av variable at . Vi vet også at tettheten til grunntilstanden er

som man kan vise direkte fra definisjonen.

Sett nå at er justert akkurat slik at , den eksakte grunntilstandstettheten. Da har vi to betingelser oppfylt:

Ved sammenlikning får vi at

Dette er en veldig banebrytende likning: Kohn og Sham viste at under antakelsen om at man kan derivere funsjonalene, så har vi en formel for det effektive potensialet! Dette understreker vi at gir en 1-elektron-Schrödingerlikning som er mye, mye lettere å løse enn det fulle -elektron-problemet.

Det er vanlig å separere ut Hartree-energien, siden denne er kjent:

der exchange-correlation-potensialet er relatert til en exchange-correlation energi,

Resultatet er en modifikasjon av Hartree-Fock-likninga som kalles Kohn-Sham-likninga, og ser ut som følger:

Sammenliknet med Hartree--Fock-likninga ser vi at leddet er erstattet med , som kalles "exchange-correlation-potensialet". Det opprinnelige leddet var kalt "exchange" og kom fra pauliprinsippet. Nå har vi en mer avansert modell, der vi har både exchange fra en mangepartikkelbølgefunksjon som ikke er en ren Slater determinant, og correlation, som ikke er tilstede i en ren Slater-determinant.

Modell-funksjonaler

Ser vi på Kohn-Sham likninga, ser vi at det som er ukjent er . Ulike modeller for dette gir opphav til ulike "DFT-funksjonaler".

Local density approximation (LDA)

I LDA antar man at XC-potensialet i hvert punkt i rommet kun avhenger av den lokale elektrontettheten , og det behandler systemet som om det var en uniform elektrongass. Det vil si at man antar at potensialet endrer seg så lite at man kan anta at elektrontettheten er så og si konstant. Dette er en tilnærming som er ganske shaky, men det gir OK resultat. Vi setter :

  • Exchange: Exchangedelen i LDA bruker den eksakte løsningen for uniform elektrongass:

    Her er en konstant som bestemmes empirisk eller analytisk, avhengig av modellen.

  • Korrelasjon: Korrelasjonsdelen i LDA er vanligvis tilpasset data fra kvante-Monte Carlo-simuleringer av uniform elektrongass. Da modelleres korrelasjonsenergien som

  • Her er en empirisk funksjon. Den funksjonalderiverte gir

LDA fungerer overraskende bra for systemer med sakte-varierende elektrontetthet, som faste stoffer eller halvledere. Molekyler og atomer er dessverre ikke slik at tettheten er saktevarierende.

Eksempel på LDA-funksjonaler:

  • Vosko-Wilk-Nusair (VWN) eller Perdew-Zunger (PZ) for korrelasjon. Dette er modeller for .

Generaliserte gradient-approksimasjoner (GGA)

Generaliserte gradient-approksimasjoner (GGA) forbedrer LDA ved å inkludere ikke bare den lokale elektrontettheten , men også gradienten . Dette gjør at funksjonalen kan ta hensyn til variasjoner i elektrontettheten på en mer effektiv måte. Intuitivt gir dette en beskrivelse av kinetisk energi.

  • Exhange: LDA modifiseres til å inkludere en avhengighet av :
    Her er en funksjon med empiriske parametre.
  • Korrelasjon: Korrelasjonsdelen korrigeres også med gradienttermer.
    Her er $g$ en funksjon.
    

GGA-funksjonaler fungerer generelt bedre for molekylære systemer enn LDA, og gir mer nøyaktige bindingslengder, energier og reaksjonsbarrierer, for eksempel.

Eksempel på GGA-funksjonaler:

  • PBE (Perdew-Burke-Ernzerhof): En av de mest populære GGA-funksjonalene. God balanse mellom nøyaktighet og beregningskostnad.
  • BLYP (Becke-Lee-Yang-Parr): Kombinerer Beckes exchange med Lee-Yang-Parr korrelasjonsfunksjonalen. Gir generelt gode geometrier, reaksjonsenergier, og vibrasjonsfrekvenser for molekyler.

Meta-Generaliserte Gradient-Approksimasjoner (meta-GGA)

Meta-GGA-funksjonaler bygger videre på GGA ved å inkludere ikke bare elektrontettheten og gradienten, men også kinetisk energitetthet. Da får vi enda mer fleksibilitet i modelleringen, men kostnadene går på den annen side opp. Den kinetiske energitettheten kommer fra Kohn-Sham-orbitalene og har uttrykket:

  • Exchange: Inkluderer en avhengighet av den kinetiske energitettheten:

der er den kinetiske energitettheten til Kohn-Sham-orbitalene.

  • Fordeler: Meta-GGA-funksjonaler er enda mer nøyaktige enn GGA-funksjonaler.
  • Ulemper: Mer beregningsmessig krevende enn GGA.

Eksempel på Meta-GGA-funksjonaler:

  • TPSS (Tao-Perdew-Staroverov-Scuseria): En av de mest brukte meta-GGA-funksjonalene.
  • SCAN (Strongly Constrained and Appropriately Normed).

Hybridfunksjonaler

Hybridfunksjonaler kombinerer noe som kalles eksakt exchange fra Hartree-Fock-teori med tilnærmet exchange og korrelasjon fra GGA- eller meta-GGA-funksjonaler. Den eksakte exchangetermen er rett og slett -leddet fra i Hartree-Fock-teori, og er inkludert for å forbedre nøyaktigheten.

  • Exchange: Man blander inn en fraksjon av Hartree--Fock exchange med en fraksjon av (meta-)GGA exchange

    der er en blandingsparameter. Her gir den deriverte av leddet som finnes i Hartree--Fock-likningene.

  • Korrelasjon: Som i GGA eller meta-GGA.

Hybridfunksjonaler er spesielt gode for organiske molekyler eller overgangsmetall-komplekser. Men, disse funksjonalene er enda mer kostbare å bruke, særlig eksakt exchange er dyrt (to-elektron-integraler.)

Eksempel på hybridfunksjonaler:

  • B3LYP (Becke, 3-parameter, Lee-Yang-Parr): En av de mest brukte hybridfunksjonalene. Beckes exchange-funksjonal er kombinert med LYP-korrelasjon.
  • PBE0: En hybridversjon av PBE-funksjonalen, med 25 % eksakt Hartree-Fock-exchange.